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Q . In [3] , we considered an approximation of a catenoid constructed from even 

truncated cones that maintains minimality in a certain sense. In this paper, 
we consider such an approximation consisting of odd truncated cones that 
maintains minimahty in the same sense. Through this procedure, we obtain 
a discrete curve approximating a catenary by exploiting the fact that it is the 
function that generates a catenoid. In this investigation, the theory of the 
Gauss hyper geometric functions plays an important role. 
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(N '■ 1 Introduction 

In this paper, we consider the problem of approximating a catenoid, which is a min- 
^ ■ imal surface, using odd truncated cones in such a manner that maintains minimality 

in a certain sense. From this surface, we also obtain an approximation of a catenary 
in the form of a polyline that retains what we regard to be the most important 
property of the catenrary. 

As a real-world illustration of the problem we consider, suppose we wish to form 
a surface from a rubber membrane. We can consider such a surface to be "stable" 
if it is difficult to deform, i.e., if it resists compression and tension. Only minimal 
surfaces are stable in this sense. For this reason, the study of minimal surfaces 
is important for industrial applications. Now, suppose that instead of a rubber 
membrane, we wish to construct a surface from plane figures. This is a common 
situation in industrial applications. A truncated cone can be constructed from 
plane figures, and thus, for such applications, it would be useful to develop methods 
for approximating surfaces of various types using truncated cones. Further, if we 
could construct surfaces that approximate minimal surfaces in such a manner that 
maintains minimality, they would be of great practical usefulness. Formulating such 
a method for constructing approximations of surfaces from truncated cones would 
also be important mathematically, because taking the limit of an infinite number 
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of truncated cones, we could obtain the parametric equations of various surfaces of 
interest. 

In this work, we approach the problem of approximating a catenoid by consid- 
ering approximations of the corresponding catenary. More precisely, we consider 
polylines forming discrete curves that converge to the catenary. Further, we choose 
these polylines in such a way that the surfaces of revolution they generate are them- 
selves minimal. Because such a polyline represents a discretization of a catenary, the 
problem we consider is also useful from a mathematical point of view. The function 
^n,2/(2m+i)(y)(0 < u < ui) (cf. Theorem 2.4) used to construct these polylines (cf. 
(2.2)) has two noteworthy characteristics: It is a rational function over the rational 
number field Q for any n and m, and it has a closed-form expression. It is interesting 
that we can approximate a catenary with polylines obtained using such functions. 

2 Main Theorem 

A catenary is the curve assumed by a hanging chain. Its form is given by the function 
Cc '■ t ccosh(t/c) for a constant c > 0. The surface of revolution generated by 
rotating a catenary Cc{t) about the t-axis is called a catenoid. Here, such a surface 
is denoted by R{Cc{t)). The catenoid has the following special property: Every non- 
planar rotationally symmetric minimal surface is congruous to a piece of a catenoid 
(cf. 3.5.1 in [2]). Here, the term "minimal" means "of mean curvature zero." Now, 
we consider catenoids -R(Cc|(_i,i)), generated by catenaries satisfying t G [—1,1], 
whose boundaries consist of two circles of radii a. If we choose a to be sufficiently 
large, then there are two such catenoids, generated by catenaries whose values of c 
we write c^, where we choose c~ < c^. It is known that the area of i?(C^+ |(_i^i)) 
is invariant with respect to infinitesimal perturbations of the surface that keep the 
boundaries fixed. The same holds for i?(C^- |(_i^i)) (cf. Theorem 1 in 2.1 of [2]). 
Moreover, it is known that the area of -R(C^+ |(_i^i)) is minimal among the set of 
surfaces possessing the same boundaries, while that of i?(C^- |(_i^i)) is not. Below, 
we consider discretizations of i?(C^± |(_i^i)). 

We call a cone whose apex is cut off by a plane parallel to its base a truncated 
cone. For xq, xi > and i > 0, let -Di/(xo, xi) be the truncated cone whose circular 
edges have radii Xq and Xi and whose height is i. Here, we do not consider the regions 
interior to the two circles of radii Xq and xi to be part of Di £{xq,Xi). Defining 

Si^e{xo,Xi) := (a;o + Xi)v^(xi -xo)^ + P, 

we see that the area of Di^i{xo, xi) is equal to n ■ S'i,^(xo, xi). For xq, xi,X2 > and 
i > 0, let D2/{xo,xi,X2) be the figure consisting of the union of Di^i{xo,xi) and 
Di^i{xi, X2) attached along the circle of radius xi. We define D 
for n > 3 similarly. The surface Dn/{xo,xi, ■ ■ ■ ,x„) is called a piecewise truncated 
conical surface with length (n; i), or simply a PTC surface with L-(n; i). Note that 
the boundary of D„^^(xo,Xi, ■ ■ ■ ,x„) consists of two circles of radii Xq and x„, and 
its area is given by 

n 
i=l 
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Next, we define 

n 

Sn/{Xo, Xi, - ■■ Xn) := ^ 5'i/(Xi_i, Xj). 

i=l 

For arbitrary fixed a,b > and n G N := {0, 1, ■ ■ ■ }, Z^.„+2/(a, Xo,Xi, - ■ ■ , b) 
is called a PTC surface with boundary condition (a, 6) and length (n + 2,i), or 
simply BCL-(a, b]n + 2; i). A PTC surface _D„+2/(a, Xq"'', x^'^'', ■ ■ ■ , xi"'', b) with BCL- 
(a, 6; n + 2; i) is said to be minimal if (xq'^'' , x^^'' , ■ ■ ■ , xi'^'') is a critical point of the 
function 

(xo,xi, • • • ,x„) H- S'„+2,£(a,xo,xi, • • •x„,6). 

Moreover, a PTC minimal surface -D„_|_2/(a, Xo°\ x^^'', ■ ■ ■ , Xn \ b) with BCL-(a, 6; 72 + 
2; is said to be stable if and only if the Hessian matrix of the above function is 
positive definite at (xq"'' , x^'^'' , ■ ■ ■ ,Xn ^). 

Studying a symmetric PTC surface D2m,i{cL, Xm-i,Xm-2, ■ ■ ■ ,Xo, - ■ ■ , Xm_2, x^-i, 
a) with BCL-(a, a; 2m; i), the following questions come to mind: 

Question 2.1. 

1. Where does this surface become minimal? 

2. Where does this surface become stable? 

We answered these questions in [3]. The answer to the first is the following: 
Answer 2.2. (cf. Theorem 1 in [3]) 

For n G N and £ > 0, there exist an explicit function gn/ix) on R>o and rjni such 
that for any m G N, the following holds: 

(1) If a > ?7m/; then the equation gmi{x) — a = has two positive solutions x^^^ 

w^th x;„ , < x+ 

(2) £'2m/(a, Xm-l,Xm-2, ■■■ ,Xq,--- , Xm-2, Xm-1, «) With X„ = gn,e{^tni,i) f'^'^ ^ = 

0, ■ ■ ■ , m — 1 are PTC minimal surfaces. 
Moreover, we have 




where Tn is a Chebyshev polynomial of the first kind, and F is the Gauss hypergeo- 
metric series. 

The answer to the second question is the following: 

Answer 2.3. (cf. Theorem 2 in [3]) Under the same conditions as in Answer 

2.2, £'2m/(a, Xm-l, Xm-2, " " " , 3^0, " " " , Xm-2, Xm-1, O) With X„ = gn/{Xa,m/) M ^ = 

0, ■ ■ ■ , m — 1 is stable. 
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In this paper, we consider the same questions for another symmetric PTC surface, 
namely D2m+i/(a, Um-i, ym-2, ■■■ ,yo,yo,--- , ym-2, Z/m-i, a) with BCL-(a, a; 2m + 
1;£). A calculation provides the following answer to the first question: 

Theorem 2.4. For n G N and £ > 0, there exist an explicit function hn/{y) on 
{y G R>o; y > ^/2} and such that for any m G N, the following holds: 



(1) If a > Um,i, then the equation hm,e{y) — a = has two positive solutions y^ 
with y-^j < y+j. 



± 



(2) D2m+l,M^ym-l,ym-2,- ■ ■ ,yO,yO,- ■ ■ ,ym-2,ym-l,a) With ?/„ = hn,i{y%^^() for 

n = 0, ■ ■ ■ , m — 1 are PTC minimal surfaces. 

Moreover, we have 

hnAy) = yVn{l + 7^ ] =yF{n + l, 



_ {2y + £)2"+^ + {2y - £)2"+^ 

where Vn is a Chebyshev polynomial of the third kind. 

The answer to the second question is 

Theorem 2.5. Under the same conditions as in Theorem 2.4, D2m+i/{0'i Vm-ii ym-2, 
■ ■ ■ , 2/0, yo, ■ ■ ■ , ym-2, ym-i, a) with yn = K/{yXm/) forn = 0,--- , m - 1 is stable. 

Now, the PTC minimal surfaces considered in Answer 2.2 with i = 1/m are 
denoted by Ei^ami/m)^ those considered in Theorem 2.4 with i = 2/(2m + 1) 
are denoted by E{yam2/{2m+i))- ^^om the fact that the areas of -R(C^+ and 
i?(C^- |(_i^i)) are both invariant with respect to infinitesimal perturbation. Answers 
2.2 and 2.3 and Theorems 2.4 and 2.5, we find that as m — )■ oo, E(x^ -, , ) and 

' ' \ a,m,l/m' 

^(<m,2/(2m+i)) ^end to while E^x^^ ^J and ^(y;;^,2/{2m+i)) tend to 

-R(C^- Therefore, we can regard the polylines whose vertices are specified by 
the sequence 



(-1' «)' ( 9m-l,l, )).■■■. (O, X^^^ ) , 



m — 1 



m 



(2.1) 



and 



2/" 



2m + l '^'™'2;4t/' 'V2m+r ""^'s;^ V^'^.^'a;^ 



(2.2) 



as good approximations of the catenaries on [—1, 1] with boundary a. 

We prove Theorem 2.4 in §3 and Theorem 2.5 in §4. In §5, by investigating 
specific catenaries, we determine the precision of the approximations of the catenary 
provided by the above polylines. 
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Remark 2.6. We showed that E{x^^-^^^) and -^(2/^.^2/(201+1))' respectively, tend 
to a catenoid by using the fact that the areas o/_R(C^±|(_i 1)) are both invariant with 
respect to infinitesimal perturbation. However, we can demonstrate the same thing 
without using this fact, since we have obtained o„ ^ix"^^ , , ) and 

771 CL^'lTh^ J. I Til ' 

^„,_2_ (^a m 2/(2m+i)) explicitly. Indeed, writing t := n/m, we have 

lim = xcosh(t/x) = Cx{t). 



Moreover, we know that x"*" -, , and x -, , tend to cj" and c„ as m — )■ 00, respec- 

' a,m,l/m a,m,l/m a a > f 

tively, because we have 

lim g^±{x) =xcosh(l/x) = Cx{l). 

In this way, by constructing PTC minimal surfaces, we can obtain the parametric 
equations of the corresponding minimal surfaces. The same is true in the case of 

3 A proof of Theorem 2.4 

We denote Si^i{s,t) by Si{s,t) for simplicity. Recall that 

Siis,t) = {s + t)^{t-sy + i^. 

Putting 

m— 1 
i=l 

The area of Z:'2.m+i/(a, Vm-i, 2/m-2, • • • , ?/o, 2/o, ■ ■ ■ , l/m-2, Vm-i, a) is expressed as 27i 
Ta,mAyo, 2/1, ■ ■ ■ , Vm-i)- Therefore, we only have to evaluate Ta^mAvo, 2/i, ■ ■ ■ , Vm-i) 
for investigating minimality and stability of this PTC surface. 
In this section, we prove Theorem 2.4. 

3.1 The case m = 1 

For a > and £ > 0, we consider the critical points of the function Ta^iAvo) (= 2/o^+ 
Si{yQ,a)). Note that 

dSi, , 2s^-2ts + e dSe, , 2t'^-2st + f 



5s v/(t-s)2 + £2' dt' ' ' ^(t-s)2+£2- 

Thus, if yo is a critical point of Ta^ij^yo), then 

^W.,,M^1^|±^.0. (3.1) 

dyo y/{a - yoY + 

Therefore, we have 

(yo - a)(%3 - Aayl + SE^yo + af) = 0. 
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Because Uq = a does not satisfy (3.1), we obtain 

« = ^o+4^. (3.2) 

If £/2 > 2/0 > in (3.2), then a is a negative number. Therefore, we regard the right- 
hand side of (3.2) as a function on {y G R>o; U > ^/2}, and denote this function by 
hi^i{y), that is, 

Afy 

Since the second derivative of hi^i{y) is 

32f y{Ay^ + 3i^) 

(%2 - £2)3 ' 

this funtion is positive and convex on {y G R>o; y > ^/2}. Moreover, because 

hm hi^e{y) = hm hi^i{y) = oo, 

j/->oo 

hi/{y) takes the unique minimal value vi/{> 0) at a point /ii,^(> ^/2). Hence, if 
a > ui^e, then there are the two solutions y^n of hi/{y) — a = 0, where y'^n > 
yUi/ > 1 £■ Therefore, ^ are the critical points of the function Ta,i,e{yo). Namely, 
^3/(a,2/Ji_£,2/Ji_£,a) are minimal. 

3.2 The case m = 2 

For a > and £ > 0, we consider the critical points of the function Taxeiuo, Hi) (= Uo^ 
+ Si{yo,yi) + Si{yi,a)), that is, we consider a point {yo,yi) satisfying 

^ (yo,2/i) = ^ (yo,yi) = 0. 

•92/0 

From the formula 

^ — (yo,yi) = 0, 
dyQ 

we see that yi = hi^i{yo). Moreover, 

. dTa,2,e , , 2yl-2ym+(^^ , 2yl-2ay, + e 

implies that 

= {2y\ - 2ym + i^f{{a - y,f + f^) _ {2yl - 2ay, + ef{{y^ - y,f + £2) 
= e{a- y^){Ayl - Aaym + 2£'yi + £'2/0 + af')- 

Here, if a = |/o, then Formula (3.3) does not hold and so we have 

[f - Ayoyi)a + Ayl + 2ey^ + ^y, = 0, 



An approximation of a catenoid 



7 



that is, 



4|/o|/i - ^2 4|/o/ii,Kl/o) - ^2 



If we put 



My) := y + 121/4^^ + 162/ (^^^ ) , (3.4) 



then /i2/(2/) is positive, convex on {y G R>o; y > ^/2} and 

hm /i2/(2/) = hm /i2/(2/) = 00. 

Thus, h2/{y) takes the unique minimal value 1^2/ at a point yU2,£- Hence, if a > i'2,e, 

then there are the two solutions ya2£ of ^2/(2/) — a = 0, where ?/a2£ > /^2,£ > 

2/~2£- Consequently, if a > U2^e, there are the two critical points {y^2i^ ^1/(2/^2^)) 

TaMvo^yi)- Namely, if a > U2/, then ^5/(0, /ii/ll/Ja/)): 1/^2/' 1/^2/' ^i/fc/))' 
are minimal. 



3.3 The case m = 3 

We consider the critical points of 

Ta,3Ayo, yi, 2/2) = yo^ + ^'^(yo, yi) + Se{yi, 2/2) + Si{y2, a) 

for a > and i > 0. If {yo,yi,y2) is a critical point of Ta^s^, then as in the case 
m = 2, we have 

yi = h^yo), 2/2 = h2,iiyo), 

and 

^ _ 4|/i + 2ey2 + £^2/1 _ 4(/i2,Kl/o))' + 2£^/i2/(2/o) + ^^/^i.Kz/o) 
42/12/2 - £2 4/ii,^(2/o)/i2/(2/o) - ^2 

_ 642/^ + 336£2?/5 + UOtyl + 7^^2/0 

- ^0 + 242/0 (^) ^- 802/0 (^) ^ + 64y. (^) ^ . 

Putting 

'^^AV) := + 24, ( + 80, (^) ' + 64, (^) ' . (3^5) 

similarly as in the case m = 2, we see that there is fi^^g > with h'^iifJ'S/) = such 
that if a > v^^g := h^^n^^, then the equation hs^y) = a has the two solutions 2/^3 ^ 
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with 2/+3^^ > /i3,^ > ya,^,t Consequently, if a > 1^3,^, then (yjg^^, /ii.KyJs,^), ^2,Kyj3/)) 
are the critical points of Ta,3Ayo,yi,y2)- 

Repeating the above argument, we see that hn,e{y) is defined as 

for n > 2 and 

hoAy) = (3-7) 
hiAy) = y + ^y^^,- (3-8) 

From Formulas (3.4), (3.5), (3.7) and (3.8), we can anticipate that hn/{y) is ex- 
pressed as, for n G N, 

KAy) = yF(n + 1, -n; I , (3.9) 

where F(q;,/3; 7; x) is called the Gauss hypergeometric series and is defined as 

^( N (tt)i(/3)i i 

and (a)i := V{a + 2)/r(Q;), and so on. Indeed, we prove this Formula (3.9) in the 
next subsection. 



3.4 Expression of hn,({y) 

In this subsection, we prove that hnAy) is expressed as (3.9) for n G N. It is obvious 
that (3.9) with n = and n = 1 are equal to (3.7) and (3.8), respectively. Therefore, 
we have only to show that Formula (3.9) satisfies (3.6) when n > 2. 
The Jacobi polynomial Pn°''^\x) is expressed as 



1 — X 

a + /3 + 1 + n, —n; a + 1; — - — 



in terms of the Gauss hypergeometric series(cf. 15.4.6 in [1]). Hence, by this relation, 
we see that the right hand side of (3.9) is equal to 

A+^^L_y (3.10) 



(1/2)/ " V V-£2y • 

In addition, by using the third kind Chebyshev polynomial VAx) which is expressed 



as 



(cf. 1.2.3 in [4]), we see that (3.10) can be expressed as 

2f 
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Consequently, we found that the right hand side of (3.9) is rewritten in terms of the 
third kind Chebyshev polynomial and denote this by hn,i{y), that is, 

KM ■=yyn{i + , , , J . (3.11) 



4r- 

So, it suffices to prove that hn/{y) satisfies the same formula as (3.6) when n > 2. 

We remark that the third kind Chebyshev polynomial satisfies the following 
relation: 

V;(x) = 2xK-i(x) - K-2(a;), (3.12) 

where Vo{x) = l,Vi{x) = 2x - l(cf. 1.2.3(1. 12a) in [4]). Then, we see that 

V^_,ix) - Vnix)Vn-2ix) = 2 - 2x (3.13) 

for n > 2. Formula (3.13) is showed in the same way as Lemma 1 in [3]. 
Rearranging (3.6), we see that the formula we should show is 

"^hn-l/iy) (hl-l/y) - hn/{y)hn-2Ay)) + (hn/iy) + 2/i„-l/(?/) + hn^2,£{y)) = 0- 

(3.14) 

Substituting (3.11) for each hi^y) in the left hand side of (3.14), and using (3.12) 
and (3.13), we see that the left hand side of (3.14) is equal to zero. Consequently, 
we have 

( 21^ \ ( 1 

huAy) = hn/iy) =yVn\l + -—r, 7^ = n + 1, -n; -; 



Moreover, Because 



F ( n + 1, — n: -: —x^ 



Ay'^-Pj ^ I ' ' '2'4|/2-£2 



/ n n \ 271+1 / I- ^ x2n+l 

(Vl + x2 + x) +(vl + x2-a;) 



2' ; 2Vl 



X 



2 



(cf. 15.1.12), we obtain 

( 1 -e \ (2y + £)2"+i + (2y-f)2"+i 
hnM = yF[n + 1, -n- -; J = 4(4^2 _ ' 

3.5 A proof of Theorem 2.4 

In this subsection, we prove Theorem 2.4. Because 



hnM = yP[- + 1, -n; -; j = , g ^,/2).(l). W - 

_^ {n + l)i fn\ y 
^ (1/2), \^) {Ay^-Py 



(3.15) 
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we have 

dy' (1/2). \ij dyH4y^-i^y 

^ (n + 1), fn\ ,, 8zym-A)y^ + 3f) 

^ (1/2), \zj (%2-£2).+2 ■ 

Therefore, hn/{y) is positive and convex on {y G R>o;y > ^/2}, and 

hm hn,i{y) = hm hn,e{y) = oo. 

Hence, there is a unique zero point fin/ of h'^iiy)- Moreover, if we put Pn/ '■ = 
hnAf^n/), then u^/ is the minimum of K,i{y). 

The role of and minimality of -D2m+i/(a, ^m-i, ym-2, ■ ■ ■ , 2/o, ^o, ■ ■ ■ , 2/m-2, 
l/m-i, a) with yn = hn/iy'^m £) for n = 0, ■ ■ ■ , m — 1 are obtained similarly as in the 
case m = 1, 2, 3. 

Remark 3.1. We saw that hn/{y) is positive and convex on {y G R>o; 1/ > ■^/2} in 
the above. We consider other properties ofhn,e{y)- 

From (3.15), we have hn+i/{y) > hn/{y) forn G N. In addition, because 

dhn.e , . in + l)i /n\ ^i d y 

dy ^ (1/2), \tj dy{Ay^-Py 

_^ {n + iy (n\ M-Sz)y^-e 
^ (1/2), \i) (V-£2)m ' 

we have h'neiv) > ^'n+iiiv) forn G N. Therefore, we obtain 

Lemma 3.2. fin+i,£ > fJ'n,e arid I'n+i/ > ^n,e forn G N. In particular, if h'^^y) > 
at a point y, then h'^ ^{y) > for n = 0,1, ■ ■ ■ ,m — 1. 

4 A proof of Theorem 2.5 

In this section, we prove Theorem 2.5. For this, we investigate whether D2m+i,^(fl, 
ym-i, ym-2, ■ ■ ■ ,yo,yo, - ■ ■ , ym-2, Vm-i, a) with ?/„ = /in/(yj„,£) fom = 0, ■ ■ ■ , m - 1 
are stable. 

We denote Si{s,t) by S{s,t) for simplicity. 
4.1 Elements of the Hessian matrix 

We give an expression of each element of the Hessian matrix of the function Ta^m/ivo, 
2/1, ■■ ■ ,2/m-i) at {y, hi^y), ■■■ , hm-i,e{y)) with a = hm,e{y). We call this Hessian 
matrix Hm{y) = (i?i,i(|/))i,j=i,2, -,m- Recalling that 



Ta,mAyo^yi^ ■ ■ ■ ^Vm-i) = 2/0^ + S'(?/o, Z/i) + S{yi,y2) H h S{ym-i,a), 
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we see that the elements of Hm{y) are expressed as 
Hi,i{y) = -^{ho,i{y),hi^e{y)), 

Hi,i{y) = -^{hi-2Ay)^hi^iAy)) + -Q^{hi~iAy)^Keiy)) 

for i = 2, 3, ■ ■ ■ , m — 1, 

Hm,m{y) = -^{hm-2,l{y)ihm-lAy)) + ^ , o) , 

Hi,i+i{y) = Hi+iAy) = Q^ihi-iAy)^Kiiy)) 



for i = 1, 2, ■ ■ ■ , m — 1, and 

if > 2. 

Now, we consider the following function: 

fmAyo,yi, ■ ■ ■ ,ym-i,ym) = yo( + S{yo,yi) + S{yi,y2) H h S{ym-i,ym) 

and define Hm{y) = {HiAy))i,j=i,2,-,m as 

HiAy) '■= ^ — ^ — hiAy), ■■■ , Kn-iAy)^ h^Ay)) 

for i,j = 1,2, ■ ■ - m. Then, we see Hi j{y) = Hij{y) except for = {m,m) and 



Hm,m{y) = -^{hm-2Ay)^hm-lAy)) + ^ (^m-l,^^) > VKZ/) ) ■ 



Moreover, we notice that liy = ?/+^_^ oi y = y',^^^, then H„^^^{y) = Hm,m{y)- There- 
fore, for investing whether I)2m+i/(a, ■ ■ ■ , l/o, 2/o, ■ ■ ■ , Z/rn-2, Z/m.-i, «) with 
y-n = hnAyamd fo^' = 0, ■ ■ ■ , m — 1 are stable, we evaluate Hm{y) instead of 

Hmiy)- 

Remark 4.1. From the construction of hm/{y), we see 

-K^iy, hAy)^ ■ ■ ■ hm-iAy)^ hmAy)) = o (4.i) 
oyi 

for < i < m — 1 . 

4.2 An expression of the determinant of Hm{y) 

We consider the determinant of Hm{y) for m > 1. Note that each det Hm{y) satisfies 
the following recurrence relation: 

det#m(i/) = ^„,„(i/)det^„„i(t/) - deiHm-2{.y) (4.2) 

for m > 3. In this subsection, we prove 
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Lemma 4.2. 

det Hmiy) = {-irHiMH2M ■ ■■H^,m+i{y)h'mAy) 

for m > 1. We remark that Hmiy) does not have Hm^m+iiy) its element. Here, 
we take Hm,m+i{y) as 



""^^'^ {y, hi,i{y), ■ ■ ■ , hm-iAy)^ hm/iy), hm+i/iy))- 



dym.-iy 

We prove this lemma by induction. 

First, we consider in the case that m = 1. Because 



dyo 

from (4.1), we have 



Observing that 

jr^{y,hiAy)) = jr^{y,hiAy)^h2Ay)) = HiAy), 
oyoyi oyoyi 

we obtain Lemma 4.2 in the case that m = 1. 

Second, we consider in the case that m = 2. Since 

hAy)^ hAy)) = O' 

oyi 

we have 

^ ^ (^^^y^ h.MMAy)) j = H,Ay) + H^Ay) ■ KAv) + ^^Ay) ■ KAv)- 

and consequently 

(-i)2#i,2(y)^2,3(y) ■ K^,{y) = - - HiAy)H2Ay) ■ KAv) 



- (^H^,2iy)y + H2,2{y)detH,{y) = det H^iy) 



from the case m = 1. 

Finally, we assume that this lemma holds for 1, 2, ■ ■ ■ m — 1. We notice that 

= T" I ^-^(y^ hiAy)^ ■ ■ ■ hm-iAy)^ h^Ay)) 

dy \oym~i 

= Hm-l,m{y) ■ Kn-2Ay) + H,n,m{y) " Kn-lAy) + ^m,m+l(y) " KnAy)' 
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Therefore, by assumption and (4.2), we obtain 

+ i-ir~'H,,2iy) ■ ■■H^^i,miy)Hm^i,miy) " Kn-2Ay) 

= Hm,m{y) det Hm~l{y) - (jlm-l,miy)j det Hm-2{y) = det Hm{y)- 

Thus, we could show Lemma 4.2. 

4.3 A proof of Theorem 2.5 

In this subsection, we prove Theorem 2.5. 
Recall that S{s, t) = {s + t)^{t- s)^ + & and 

dsdt^'' ^ ((t_,)2 + £2)3/2' 

Hi,i+i{y) = Hi^i+i{y) = Q^{hi-i,e{y),hi^i{y)) 

for i = I,-- - ,m. Therefore, every Hii^i{y) is negative on {y G R>o;y > ^/2} 
and thus, the sign of det Hm{y) is equal to the sign of h'^^{y) by Lemma 4.2. The 
following lemma is well known: 

Lemma 4.3. A symmetric m x m matrix A = (Aj)i.j=i,2,- -,m is positive definite if 
and only if det A^^^ > for any n = 1, 2, ■ ■ ■ , m, where A^*^^ := {Aij)ij=i^2,- ,n- 

Note that Hm\y) is equal to Hn{y)- If 2/ > A^m/, i-e., y = ya,m,e^ then det Hm{y) = 
det Hm{y) > and detHn{y) = det Hn{y) > for n = 1,2, ■■■ ,m — 1 from 
Lemma 3.2. Hence, from Lemma 4.3, we see that Hmiyami) positive definite 
and -D2m+i,K«5 2/m-i,2/m-2, ■■■ ,yo,yo,--- ,ym-2,ym-i,a) with yn = K/iy^^^A for 
n = 0, ■ ■ ■ , m — 1 are stable for a > I'm/- 

5 Approximations of catenaries 

In this section, we see that the polylines (2.2) considerably approximate catenaries, 
by calculating some cases numerically. 

Recall a catenary Cc(t) := ccosh(t/c). The function c i— >■ ccosh(l/c) is positive, 
convex on R>o and takes the unique minimum rjoo '■= 1.5088 ■ ■ ■ at ^oo 0.8335 ■ ■ ■ . 
Thus, if a > r/oo, then there are two positive numbers with c~ < ^oo < such 
that c^cosh(l/c^) = a. 

Before approximating catenaries, we investigate a relation between t]^ and 
'^m,2/{2m.+i) in the following subsection. 
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5.1 A relation between r/oo and ^'^,2/(2™+!) 
In this subsection, we show the following lemma: 
Lemma 5.1. r]oo > i^m,2/(2m+i) for m>l. 
It suffices to prove that 

ycosh{l/y) - hm,2/{2m+i){y) > (5.1) 

for y > I/V2 because the first derivative of ycosh(l/y) at y = l/-\/2 is less than 
zero. The left hand side of the Formula (5.1) is equal to 

?/cosh ( — ) — yF ( m + 1, — m; 



(y) V"" ' ^' 2' (2m + 1)V _ 1 

t^(l)2.y2-^ (1/2). ^zy'((2m + l)V-l)^ 

E°° / 1 (m + l)j /m\ y \ 

1(1)2.1/2-1 (1/2), ^zy'((2m + l)V-l)V 



i=0 

00 



=: 2^a(m,2,?/). 



i=0 



Note that a(m, 0, y) = 0. Now, we investigate a(m, i, y). First, we evaluate a(m, 1, y) 
+a{m,2,y). This is equal to 

8m(m + l)y^{2y^ - l)i3y^ + 1) + {y - l^iy + l^iUy^ + 1) 
24?/3((2m + 1)V - 1)^ ' 

Therefore, if ?/ > l/"\/2, then a(m, 1,?/) + a{m,2,y) is positive. Next, we evaluate 
a{m,i,y) for z > 3. Because 

1 1 

> 



(2m+l)2?/2_i 4(m + 2)(m- l)?/2 
for ?/ > 1/3, 

1 (m + l)j / m\ y 



(1)2.1/2*"^ (l/2)i \i J A'{m + 2y{m-iyy^' 
^_ im + i)\ 



a{m,i,y) > 



(1)2. ?/2^-^ V (m - i)!(m + 2)*(m - 1)* 

for ?/ > 1 / -\/2. We show 

1 (m + z)! 

1 > -77 — =: b(m, i) 

{m-i)\{m + 2y{m-iy ^ ' 

for z > 3. We easily find that 6(m, i) > b{m, i + 1) for z > 3. Therefore, we have 
only to show 1 > b{m, 3). This is shown by 

b[m, 3) — 



(m-3)!(m + 2)3(m- 1)^ 
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_ m(m - 2) ■ (m + l)(m + 3) (m - 1)^ ■ (m + 2)^ _ ^ 
~ (m + 2)2(m- 1)2 ^ (m + 2)2(m-l)2 

Hence, a{m,i,y) is positive for z > 3 and y > l/-\/2. Collecting the above results, 
we see that (5.1) is correct for y > l/\/2. 

5.2 Approximations of catenaries 

By Lemma 5.1, if a > ?7oo, there are ya,m,2/(2m+i)± for m > 1, that is, we can 
approximate C^±{t) by the polylines (2.2). 
For example, if a = 2, then 

yti,2/3 ^ 1-7338, 

^3,2/7 ~ 1-7035, /ll,2/7(<3,2/7) ^ 1-7518, /i2,2/7(<3,2/7) ~ 1-8497, ■ ■ ■ 

C2 ~ 0.4701 and 

C^+(l/3) ^ 1.7295, 

C^+(l/5) ^ 1.7084, C,+ (3/5) ^ 1.8038, 

^4(1/7) ^ 1.7026, C,+ (3/7) ^ 1.7510, ^4(5/7) ^ 1.8492, ■ ■ ■ . 

Moreover, 

?/2~l,2/3 ~ 0.5150, 

l/-2,2/5 ^ 0.4856, /ii,2/5(2/2:2,2/5) ^ 0.8823, 

^2~3,2/7 ~ 0.4779, /il,2/7(Z/2~3,2/7) ~ 0.6655, /i2,2/7(y2~3,2/7) ~ 1-1141, " " " 

^ 1.6966 and 

C,-(l/3)^ 0.5933, 

C^-(l/5) ^ 0.5133, C^-(3/5) ^ 0.9078, 

C,-(l/7) ^ 0.4920, C,-(3/7) ^ 0.6794, C,-(5/7) ^ 1.1254, ■ ■ ■ . 

So, we see that the polylines (2.2) form nearly equal catenaries C*c*|{ 1 1) large 
enough. We draw two polylines whose vertices are specified by the sequence 



2m + l V2m+1 2^^1+1 V^«.™.ra 

in cases that (a, m) = (2,5) and (2, 10) below: 
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